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I. INTRODUCTION 



Recently phenomena involving interactions of hadrons in small-size configurations have 
been intensively discussed both in relation with the phenomenon of color transparency and 
vector meson electroproduction observed at HERA energies. There is also a deep relation 
between presence of the weakly interacting small size configurations in hadrons and phe- 
nomenon of cross section fluctuations in the interactions of the hadrons which manifests 
itself in the inelastic coherent diffraction processes: h + N(A) —>■ X + N(A), see Ref. [|TJ. 
In this paper, we focus on the systematic derivation of the formulae for the interaction of 
the color singlet qq pair having a small transverse size with a hadron target. Then, we use 
these formulae to calculate the probability of the distribution for the interaction of a photon 
and a pion with a target for small interaction cross sections. Although some of equations 
deduced in the paper existed before no derivations with analysis of their accuracy has been 
presented. 

The paper is organized as follows. In Section 2, we consider the virtual forward Compton 
amplitude in the small- x region where it is dominated by the photon-gluon scattering sub- 
process. We outline there a derivation of the basic formula expressing the total cross section 
(X-y^T 1 as a convolution of the gluon distribution amplitude G T {x,Q 2 ) and the jg scattering 
cross section. In Section 3, we write down the 7*7 cross section in terms of the qq light-cone 
wave functions of the virtual photons. In the next section, we calculate the cross section 
distribution P 7 »(cr) for the virtual photon. In Section 5, we discuss the quark-hadron dual- 
ity interplay between the perturbative free-quark results and contributions due to low-lying 
resonances. Finally, in Section 6, we calculate the cross section distribution for the pion 
P 7r (cx) in the small cross section limit where it is governed by qq configurations having small 
spatial size. Basing on QCD evolution equation we evaluate also the functional dependence 
of P n (cr) — > on a and on the incident energy. 
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II. HARD 7*T TOTAL CROSS SECTION AND THE INTERACTION OF SMALL 

SIZE CONFIGURATIONS. 



Let us consider a particular contribution into the 7*T cross section corresponding to a 
transformation of the virtual photon when 7* converts into a QQ pair with quarks having a 
large relative transverse momentum. Usually, this contribution is written as a convolution 
of the infinite momentum frame wave function of the target with the pQCD calculable 
coefficient function describing the short- distance propagation of the particles between two 
virtual photon vertices. Our aim is to express the relevant coefficient function in terms of the 
light-cone wave functions of the virtual photons as viewed from the reference frame where 
the target is at rest. The contribution we are interested in is given by the sum of diagrams 
shown in Fig.|l|. 

The lower blob corresponds to the gluon distribution in the target. It is convenient to 
parameterize the gluon momentum k in terms of the Sudakov variables 

k = -aq' + Pp' + kt d 4 k = ^dad(3d 2 k t . (2.1) 

Here q' and p' are light-like momenta related to p, q by 

q = q ' + 2W) p ' P = p ' + W7) q '-' (2 ' 2) 

2{„) = 2(pV) + ■ (2.3) 

For our goals, the most interesting region is that of small values of the Bjorken parameter: 
x = —q 2 /2(pq) — > 0, where we may safely approximate 2(p'q') = s. 

In the d A k integral, the region k 2 ~ Q 2 corresponds to the next-order a s correction, so we 
will take into account only the contribution of the region k\ <^Q 2 . This corresponds to the 
leading a s logQ 2 /A 2 approximation in which the 0(a s ) corrections are neglected. In this 
kinematical region, the contribution of the diagrams shown in Fig. 1 can be considerably 
simplified. The essential region of integration is 



I k 2 1 = 1 a(3s + k 2 t |<Q 2 . (2.4) 

Note that [3s ~ Q 2 oc (mass) 2 of the qq state produced by 7*. Hence, 

a < 1 (2.5) 

is the essential region of integration over a. It is convenient to write the propagator d^(k) / k 2 
of the exchanged gluon in the light-cone gauge q'^A^ = 0, in which 



(W) 

It can be shown (cf. ||) that in our case the k^q'^ part of the propagator dominates. Using 
eq. Q2.5p , the dominant part of the gluon propagator can be further simplified: 

tLa(jfe) ~ (2.6) 

In other words, it is sufficient to take into account only the longitudinal polarization 
of the exchanged gluons. Indeed, let us estimate the contribution due to exchange of a 
transversely polarized gluon: 

Here, Tff£) is the imaginary part of the amplitude of the 7* scattering off a gluon given by 
the lowest-order Feynman diagrams and T^ ± \ ± is that for the gluon scattering off a target T. 
Using the fact that, at high energies, the Feynman amplitude of processes due to exchange 
by two elementary fermions tends to constant || we obtain: 

((3s) 2 



J Tf ±x± d({3s) cx ^ ( g 2) . 



In this estimate we use also scaling over Q 2 in the box diagram. The amplitude due to the 
exchange by two vector particles increases like s, and we have 



J T^ x± das oc {as) 2+n . 



Here n > 0, since, according to the QCD evolution equations, the deep inelastic amplitudes 
increase with energy in region of applicability of the perturbative QCD. 



As a result of this power counting estimate we obtain: 

f (f3s) 2 (as) 2+n d 2 k t n ,~ , n 

fl(fc t 2 < ^)(as) n oc (as) n . (2.8) 



Here, we substituted a(3s ~ k t . Thus, due to the presence of the factor a in eq. (|2.8|) , in the 
leading a s logQ 2 / X 2 approximation, the contribution due to the exchange of a transversely 
polarized gluon is negligible compared to the contribution of the longitudinal polarization 
specified by eq.( |2.6| ). We use here the observation that in QCD the power n characteriz- 
ing the energy dependence of the amplitude is the same for scattering of transversely and 
longitudinally polarized gluons. 

Using the gluon propagator in the form given by eq. ( |2.6|) , we get the following expression 
for the total contribution of the diagrams shown in Fig. 1: 

ImM = / J^H Wf 2 ' ^x {P) ^^ T) d^k)d xl {k) (2.9) 

Here T°^ P ' = T^ x (^*g — > QQ) is the sum of the box diagrams describing the *y*g scattering 
and T^j T ^ is the amplitude of the gluon scattering off the target T. 

Using the dominance of the longitudinal gluon polarization ( |2.6|) and incorporating 



eq. (|2.1|) we can rewrite eq. (|2.9|) as 



ImM_ f sdad(3d 2 k t 2ImTg (P) p M p A 4Imrg (T) q m 



2(27r) 4 (fc 2 ) 2 4(pqy 
Now we will use the fact that 



(2.10) 



T$ P K = T^ P h x = 0, (2.11) 

since the box diagram contains no gluons and, therefore, the Ward identities in this approxi- 
mation are the same as in an Abelian gauge theory (a, b are the color indices). From eq. Q2.1| ) 
and (|2.11 ) it follows that 

lmT^ b x {P) p^ Px _ ImT^fc^ 



4(pq) 2 (f3s) 2 (2 ' 12) 



and 

ImM - /•rf.g JSffifl / 4 ^ ImT g (T W 

s ~y ^ (/3s) 2 J (2tt) 4 (A; 2 ) 2 S 2 • 1 ' 6) 

It is useful to define the cross section of 7* scattering off a gluon g averaged over the gluon 
color: 

6a b -s'-a(Y9^qq) = l £ > (2.14) 

where s' is the invariant mass of the produced qq system: s' — (k + q) 2 ~ j3s — Q 2 . Thus, 
ImM fd(3 _ f sdad 2 k t 4ImT_^ >qp,q- x 



a ^ T =—=j j a ^ 9 -> ^ y (2^)W^ ? — ^ — • (2 - 15) 

In the leading a s logQ 2 approximation, we can substitute k 2 by k 2 . Comparing our result 
with the QCD-improved parton model expression for the production of heavy quarks (see 
e.g., ||), we observe that 

/ f^E4Im^ (T> te« = PG T (fi,Q\ (2.16) 

where Gt is the gluon distribution in a target T. This gives 

a rT = Ja rg ^j-{f3G T (P,Q 2 )}. (2.17) 

The first argument of G{(3, Q 2 ) is (3 — ® + s M . Here M is the mass of the produced qq pair, 
which is typically of the order of Q. Hence, the essential region of integration is /3 ~ x. As 
the evolution scale for the gluon distribution function, we take Q 2 . Of course, higher-order 
a s corrections may change Q 2 by some numerical factor. This scale-fixing ambiguity is a 
usual feature of the leading a s log (^-approximation. 



III. LIGHT-CONE WAVE FUNCTIONS AND oy G 

Now let us express <r 7 * 9 in terms of the light-cone wave functions of the virtual photon. 
To this end, we write down the four-momenta ri(r2) of quark (antiquark) in the box in 
terms of the light-cone variables 7*1 = {r^rf ,r t } with rf = rjq + and take the integral over 
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r l by residue. Introducing the lowest-order perturbative qq light-cone wave functions of the 
virtual photon 

i (31) 

we obtain the following expression for the sum of the box diagrams: 



•ab(p) 

das 



ImT /t> P^Px 



s 2 



e2 9s / /o r t){2^(v, r *) - ^7* fa r < + fc t) - Vvfa r * - fct)}^^, (3.2) 



di]d 2 r t 

where g 2 is the QCD coupling constant and F a = 4f, A a being the Gell-Mann matrices of 
the SU(3) C group in the fundamental representation. 

It is convenient to rewrite this formula in the impact parameter space: 

Mz,r t ) = J M^b)e irtb d 2 b. (3.3) 

Then 

/ ^ ^pqy^ = I ^(-^)^ s 2 {-[2-e^- e -^]Tr(F a F fe )}. (3.4) 



Within the leading a s logQ 2 approximation, to obtain eq.( |2.12|) , it is necessary to decompose 



exponent into a power series over (ktb) and to keep terms up to the second order in k\ 
Combining eqs. (|3.3|), (|2.10|) , and ( |2.16| ), we obtain 



^T = e 2 /^(r/,5)^iV c |^>^Tr^}G T (x,A/6 2 ). (3.5) 

Here, factor A can be estimated from analysis of (Tl{i*N) cross section. Since the gluon 
density increases when x decreases, A slowly increases with decrease of x ||. For x ~ 1CT 3 , 
A ~ 9. 

It is instructive to represent ct 7 *t in the form 

a rT = e*J^b)^N c .a?(b>). (3.6) 

Here, of 9 is the cross section for the interaction of a colorless small transverse size gg-pair 
with the target T: 



of = dl^^Tr ^pj xG T (x, A/6 2 ). (3.7) 

This expression was obtained originally in As usual, N c is the number of colors, and 

the Casimir operator of the SU(3) group in the fundamental representation can be easily 
calculated: 

-—TrF 2 = -TrFn = -. (3.8) 
8N C 3 3 6 ( J 

Combining all the numbers together, we finally obtain: 



7T 2 

< = -b 2 

T 3 



xG T (x,X/b 2 ) a s (X/b 2 ). (3.9) 



Here b = (b q — bq). This formula describes the essence of the color transparency (CT) 
phenomenon (cf. discussion in 0): qq configuration of a small spatial size has a small 
interaction cross section. However, for sufficiently small x, the interaction becomes strong 
due to the formation of the soft gluon field. In this respect, eq. ( |3.5| ) predicts the interaction 
of a small size configuration which is qualitatively different from that of the models of 
F. Low [|nj and J. Gunion and D.Soper |TT[|. The fact that erf? is proportional to the 



gluon distribution in eq. (|3.7| ) increasing in the small- x region, has important experimental 
consequences, e.g. it makes it possible to observe the small-size quark configurations at 
HERA in the electroproduction of vector mesons at small x. In fact, eq. ( |3.5|) can be inferred 
from a formula derived in ref. within a model approximation to QCD. Using some simple 
tricks, one can also obtain eq. ( |3.5|) from a formula obtained in |13| within the leading a s In x 
approximation of QCD combined with some bold assumptions concerning the parton model 
structure. 

Using ( |3.2| ), we can calculate distribution over cross section for the fast photon or pion 
projectible for small a (c/. 0). 

IV. DISTRIBUTION OF P rN (a) FOR THE PHOTON PROJECTILE. 

In the previous section we have derived eq. (^.5[) which expresses the <t 7 *t cross section 



in terms of the light-cone wave functions of the virtual photon 7*. This formula gives us the 
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possibility to calculate another useful quantity - distribution over cross section 

P 7 * T (cr). By definition, the differential probability that the virtual photon 7* interacts 
with the target T with the cross section a. In other words, the experimentally observable 
total cross section in terms of P(cr) is given by 

<y-i*N = / -P 7 *Ar(cr)(T(i(T. (4.1) 



In refs. [|b|,[L5|], it has been suggested to represent the cross section o in terms of the 
eigenstates of the S'-matrix. In the case of small a, as a result of color screening and 
asymptotic freedom, the scattering state is a qq pair. So, the contribution of small a has 



the form of eq. ( |3.6|) . Using eq. ( |3.7| ), we can write: 



Let us define 



ea rN = e 2 / ^ (77, b) ^ N ^^~da . (4.2) 



P rN (a^0) = J e 2 ^{v,b)—N c —, (4.3) 

where ip^ij], b) is given by eqs. (|3.1|) and (p.3|) . It is implied here that the functional depen- 
dence of b on a in eq.( |4.3|) should be calculated from eq.( [3.7|) . Now, we can rewrite eq . (|4. 2|) 
in the form of eq.( |4.1|) . Though our derivation is applicable for the interactions with small 
a, eq. ( |4. 1|) has a more general nature. In fact, it has been understood long ago [p^4Vfl6| that 
many features of the interaction of a fast projectile can be described in terms of distribution 
over cross section. An important advantage of such a quantity is that it accurately takes 
into account diffractive processes. Some properties of P(cr) have been discussed in detail in 
However, for our purposes, it is sufficient to consider P 7 .jv(cr) in the limit a — ► 0. 
In general, eq. ( |4.3| ) predicts a rather involved dependence of P(cr) on In a at small a. 
However, this dependence can be easily calculated using QCD evolution equations. The 
distinctive feature of eq.(|4.3|) is that 



P 7 »at(<t — > 0) |cr^o~ — up to log(cr/cr ) terms. (4.4) 
a 
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V. TRANSITION TO MESONS 



The perturbative version of the virtual photon wave function ip^(r], r t ) (3.1) can be writ- 
ten through a dispersion integral 

*Pn{v,n) = - [ ^f{K]7],r t ) f K n2 (5.1) 
n Jo H /t z + Q z 

where 

is the wave function of a non-interacting gg-pair with invariant mass k. The interaction 
between the quarks modifies the virtual photon wave function ip^f], r t ) — > ^(77, r t ), and the 
dispersion representation 

j poo _ cIk 2 

Vv(^ r i) = -/ ^l q (^V,r t ) 2 . n2 (5-3) 

TT JO K z + Q A 



for the "exact" wave function ^^(r/, r 4 ) is in terms of the modified spectral density 
i>p (k] t\, r t ) in which, instead of the free-quark approximation ip^^n; r\, r t ), one has a sum 
over resonances, the p-meson being the dominant feature in the 1ow-k region: 

77, r t ) - ^ (k; r/, r t ) = ^fa, r t )5( K 2 - mj) + r?, r t ) (5.4) 

where g p is the magnitude of the p-state projection onto the electromagnetic current. At large 
n, the resonances are wide, and their sum rapidly approaches the free-quark value, i.e., one 
has a perfect quark-hadron duality^} For sufficiently large Q 2 , the dispersion integral ( |5.3|) 
is dominated by higher states, and the free-quark approximation is completely justified. 
Decreasing Q 2 , one would observe mismatch between the free-quark calculation and the 



Note, that since the large-K behaviour of tp l } adr (n;ri,rt) coincides with that of ipff^K) rj, rt), the 



dispersion integral in eq.( |5.3D has the same convergence properties as that in eq. fl5.1|) , i.e., no 



subtractions are needed in eq.(5.3) 
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dispersion integral over the resonances. Such a situation is well known from QCD sum 
rules: the difference between the resonance and free-quark spectra is described by power 
corrections (1/Q 2 ) N . The usual procedure is to approximate the higher states by the free- 
quark contribution ("first resonance plus continuum" model) 

^ i9herstates {K;r),r t ) = 9(k 2 > s^f^n) 

where Sq is the effective threshold for higher resonances in the p-channel and then fix its 
value by the requirement of the best agreement between the two sides of the resulting sum 
rule 

After fixing Sq from the magnitude of the power corrections A^/(Q 2 ) N , one can take the 
limit Q 2 — ► oo to get the local duality relation 

Kg P V»{v, n) = j q v, r *) dn 2 . (5.6) 

In other words, the p-meson wave function in such an approach is dual to the free-quark 
wave functions integrated over the duality interval < n 2 < Sq. 

For the forward virtual Compton amplitude, the dispersion representation can be applied 
both for the initial and "final" virtual photon. However, taking only the p-meson contri- 
bution in the dispersion integral for the final state, one naturally obtains the amplitude 
for the 7*T — > pT transition considered in ref. 0. Furthermore, picking out the p-meson 
contribution in both dispersion integrals one would get the amplitude for the pT — > pT 
scattering. This idea can be also used to study the pion diffractive electroproduction and 
the pion diffractive scattering. 

VI. CALCULATION OF P nN (a 0). 

To analyze the pion scattering, we substitute the electromagnetic current by the axial 
current in the original amplitude, i.e., simply add 75 in the current vertices. For massless 
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quarks, the final result has the same structure as that for the vector current. Of course, the 
gg-pair wave function would have an extra 75, and the vertex factor analogous to that in 
eq.(5.2) is 

U{xP + )>yyV((l-x)P + ) 

where P is the 4- momentum associated with the axial current. 

The projection of a single-pion state onto the axial current is specified by the tc — > fiu 
decay constant f n : 

(0\J^\n,P) = V2f n P,. (6.2) 

Hence, we should extract the amplitude ~ P^P P corresponding to the longitudinal polariza- 
tion of the axial current. Again, the transition from the virtual amplitude for the currents 
to that involving the pion can be understood in terms of the dispersion representation and 
quark-hadron duality. In other words, below the effective higher state threshold Sq, one 
should substitute the free-quark contribution by that due to the pion pole: 

<(«;i7,r t ) - 17, r t ) = q, (f*Mv, n)*(« 2 - ml) + 6(k 2 > a^( K] T,,r t )) . (6.3) 

The local duality prescription gives a correctly normalized wave function provided that 
Sq = 16ir 2 f 2 m 0.67 GeV 2 . Of course, one can use a pion wave function different from that 
given by the local duality. However, the duality considerations justify the use of the effective 



two-body wave function (see ||18|| ). 

The actual calculation consists of the same steps as those leading to eq. fl3~E|) . For a 
small-size configuration, we get the following contribution 5a n N into the scattering cross 
section: 

5a nN = J \^(ri,b)\ 2 ^N c a%(b 2 ). (6.4) 

Effectively, the vertex eip^y/N^ is substituted by the pion wave function. Rewriting Sct^n 
as 
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Sa nN = P nN (a)da 



dtf_ 
da 



J \^( v ,b)\ 2 ^ada , 



(6.5) 



we obtain: 



p^(^o) = ^/|^M-o)| 2 ^, 



(6.6) 



where a{b 2 ) is given by eq. (|3.9| ). 

Thus, P vr Ar(<7 — > 0) is determined by the pion wave function at the origin of the impact 
parameter space, or, what is the same, by the integral of the momentum wave function 
i>Tr(Vi r t) over a h transverse momenta r t . This integral formally gives the pion distribution 
amplitude 



However, in QCD (and in any theory with dimensionless coupling constant), this integral 
diverges. The standard procedure is to supplement the integral with some renormalization 
prescription characterized by a cut-off parameter /z, i.e., <p n (rj) — > <p n (ri,{i). In fact, the 
Fourier transformation from the momentum to the impact parameter space 



for small b can also be treated as a particular cut-off prescription with 1/6 playing the role 
of the renormalization parameter fi. In the b — > limit, one encounters the singular logfe 2 
terms. It is exactly the logarithms which generate the evolution of the pion distribution 
amplitude. Summing the logarithms by the renormalization group methods gives, for small 
b: 



where r](l — r])C^ 2 (2ri — 1) are the eigenfunctions of the evolution kernel (C^ 2 (2rj — 1) 
being the Gegenbauer polynomials), the anomalous dimensions 7„ are its eigenvalues and 
fio is the one-loop QCD /^-function coefficient. The &o-parameter characterizes the effective 
onset of the perturbative evolution. The coefficients a n are the Gegenbauer moments of the 






(6.7) 
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pion wave function at this scale. Note that the anomalous dimension of the axial current 
vanishes (70 = 0) and all other 7 n 's are positive. Hence, after the renormalization group 
improvement, the limit 6 — ^ is well-defined in this case and 

^,6 = 0) = -s/4»7r/ w 77(l-77), (6.8) 

where f w = 92 MeV. The absolute normalization of the pion wave function for b = is fixed 
by the matrix element of the axial current: 



or in the impact parameter space (see eq 

f^ b = 0)d £=7k- (610) 

In other words, for the pion, the singular log 6 terms sum into harmless (l/log6 2 A 2 ) 7n//2/3 ° 
factors vanishing in the b — > limit. As a result, the ^-dependence of the pion wave function 
ip n (r), b) in the formal b — > limit always assumes its asymptotic form ip v {r],b) ~ rj(l —77), 
irrespectively of its shape at the scale b . It is natural to expect that b is related to the scale 
characterizing the magnitude of the nonperturbative momentum distribution in the pion. 
The momentum scale fiQ = y/s^ ~ 0.8 GeV suggested by the local duality is rather large, 
and there may exist a transitional region of distances b ~ 60 small compared to the pion size 
but not small enough to produce sizable perturbative evolution effects. In this case, one can 
try the ^-dependences of ^(t], b) different from the asymptotic form. In fact, the integral 

1 = J |WM)I 2 ^ (6.11) 

is rather insensitive to the evolution effects. If we take the asymptotic wave function ( p.8| ) 
then 

Assuming that, at the scale b = bo, the 77-dependence of the pion wave function corre- 
sponds to the Chernyak-Zhitnitsky |19| ansatz 



14 



^ Z ( V , b = b ) = 5V487T/; 77(1 - V ){1 - 2 V ) 2 , (6.13) 

we obtain 

PA°,b = h) = ^fi^. (6.14) 

Thus, in this case the evolution would decrease the integral / by ~ 20% when b changes 
from b to 0. Taking the asymptotic result, we get: 

6 f 2 

P>-0) = - J f (6.15) 

5 a s x&N{x, A/cr) 

Distribution P n N(a) was determined in Ref. from the analysis of the soft diffractive 
processes for E n « 200 GeK, see solid curves in Fig. 2. In the limit a <C (a), we can compare 
this result with eq. (|6.14j ). The applicability region of this equation is restricted by several 



conditions. First, x e ff should be small enough so that the average longitudinal distances 
in the scattering process l/2rriNX are larger than the nucleon size, which corresponds to 
x < 0.05. Furthermore, the virtualities in the process should be large enough so that one 
can apply pQCD which corresponds to the requirement Qlff ~ 1 — 2GeV 2 . In our analysis 
we also neglect the 6-dependence of the wave function of the qq component at large b (this is a 
higher twist effect), which restricts consideration to b < 0.5 fm. In the numerical calculation, 



we use the GRV parameterization JSC] since it describes well the parton distributions down 
to Q 2 ~ 1.5 GeV 2 . We present results both for the leading and next-to-leading order GRV 
parameterizations, see dashed curves in Fig.2. Difference between LO and NLO results 
illustrates range of uncertainties of the current analysis. One can see that the results of our 
calculations are in qualitative agreement with the phenomenological results of |7j]. 

Another interesting feature of our results is a substantial energy dependence of P(a < 
(a)) on the incident energy due to a fast increase of xGn(x,Q 2 ) with the decrease of x, 
see Fig. 3. This reflects the fact that the probability of point-like configurations in hadrons 
decreases with the increase of energy. Further diffractive data (preferably at higher energies) 
are necessary to get better information about P^^^a). 
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Since the existence of configurations with small spatial size has been confirmed experi- 
mentally in the energy dependence and absolute value of cross section of electroproduction 
of vector mesons, we consider the above result as a reflection of soft matching between 
nonperturbative and pQCD regimes. 

VII. SUMMARY AND CONCLUSIONS 

In this paper, we applied a pQCD approach to describe the basic features of the high- 
energy interactions of a small-size qq configurations with a hadron target. This interaction 
is proportional to the gluon distribution function Gt(x,Q 2 ) of the target and, hence, the 
cross section is enhanced in the small-x region. The qq configuration can be described by the 
wave functions whose particular form is determined by the projection of the initial particle 
(7*, p or 7r) onto the qq component. For small a, we calculated the cross section distribution 
P 7r (cr) for the pion and demonstrated that it is rather insensitive to the specific form of the 
pion distribution amplitude. 
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FIG. 1. Leading small- x contribution to the forward virtual Compton amplitude. 
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FIG. 2. Comparison of P 7rp (o") calculated in pQCD using eq.( |6.14[) ) and GRV parameterizations 
1 20] of the gluon density and fits based on the analysis of the soft diffraction data Q. 
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FIG. 3. Incident momentum dependence of P np (a) for small a calculated using eq. flOU ) and 
GRV NLO parameterization [|(]] of the gluon density. 
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